Ferromagnetism in a dilute magnetic semiconductor — 
Generalized RKKY interaction and spin-wave excitations 
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Carrier-mediated ferromagnetism in a dilute magnetic 
semiconductor has been studied using i) a single-impurity 
based generalized RKKY approach which goes beyond lin- 
ear response theory, and ii) a mean-field-plus-spin-fluctuation 
(MF-I-SF) approach within a (purely fermionic) Hubbard- 
model representation of the magnetic impurities, which in- 
corporates dynamical effects associated with finite frequency 
spin correlations in the ordered state. Due to a competi- 
tion between the magnitude of the carrier spin polarization 
and its oscillation length scale, the ferromagnetic spin cou- 
pling is found to be optimized with respect to both hole dop- 
ing concentration and impurity-carrier spin coupling energy 
J (or equivalently U). The ferromagnetic transition tempera- 
ture Tc, deteremined within the spin-fiuctuation theory, cor- 
responds closely with the observed Tc values. Positional dis- 
order of magnetic impurities causes significant stiffening of 
the high-energy spin-wave modes. We also explicitly study 
the stability/instability of the mean- field ferromagnetic state, 
which highlights the role of competing AF interactions caus- 
ing spin twisting and noncollinear ferromagnetic ordering. 



I. INTRODUCTION 

The discovery of ferromagnetism in Mn-doped III- 
V semiconductors such as p-type Ini-^Mn^^As, [1] and 
Gai-ajMn^^As, [2] with a highest transition temperature 
(Tc) of UOK for Mn concentration x = 0.053, [3] has led 
to considerable interest in these dilute magnetic semicon- 
ductors (DMS). The successful search for ferromagnetic 
ordering above room temperature in Gai_xMna;N, [4,5] 
with a highest reported value of 940K, [6] has added 
a new dimension to the interest. 

Besides their potential applications in semiconductor 
devices such as optical isolators, magnetic sensors, non- 
volatile memories seamlessly integrated into semiconduc- 
tor circuits etc., and possibilities in photonics and high 
power electronics, attention has also been focussed on 
the fundamental mechanism and nature of the ferromag- 
netic state, and the possibility of studying new magnetic 
cooperative phenomena such as spin-dependent tunnel- 
ing, magnetoresistance, spin-dependent light emission 
etc. in semiconductor heterostructures arising from the 
new (spin) degrees of freedom. 

The double-exchange model, involving the interaction 
— JSi.di between the magnetic impurity spin Si and the 
electron spin cji, has been the starting point in nearly 



all theoretical studies, and we first review the emerging 
physical picture and the different approaches employed. 

Long range ferromagnetic interaction between the S — 
5/2 Mn++ ions is mediated, in the mean-field (Zener 
model) picture, [7-13] by a uniform itinerant-carrier spin 
polarization, which is caused, in turn, by an effective uni- 
form magnetic field, resulting from site-averaging (virtual 
crystal approximation) of the local impurity fields. In 
the weak- field limit [xJS « ep), the carrier spin po- 
larization is proportional to the Pauli susceptibility xp; 
and the transition temperature (Tc ~ xJ^xp) is therefore 
proportional to the Mn concentration x, J^, the carrier 
effective mass m*, and iV(eF) ~ P^^^^ where p is the 
hole concentration. In DMS, p is a only a small fraction 
(/) of X due to large compensation by As antisite de- 
fects. Therefore, the Fermi energy ep ~ Wp^^^ itself is 
quite small compared to the bandwidth W, and hence the 
weak-field limit is vahd only for x « {W/ JS)^f'^. Va- 
lence band spin splitting comparable in size to the Fermi 
energy has been confirmed experimentally. [14] 

Dynamical correlations in the ordered state have been 
studied within a path-integral formulation in which the 
itinerant carriers are integrated out and the effective ac- 
tion for the impurity spins is expanded up to quadratic 
order (non-interacting spin- wave approximation). [15] In 
contrast to the MF results, the spin stiffness (and hence 
Tc) is independent of J and inversely proportional to 
m. Other approaches incorporating dynamics include the 
dynamical mean field theory, [16,17] in which the local 
charge and spin fluctuations are included but long-range 
spin-wave excitations are neglected, and a RPA-level 
spin-fluctuation approach in which Mn disorder is treated 
within the coherent potential approximation (CPA). [18] 

While the positional disorder of Mn ions is not taken 
into account in the virtual crystal approximation (VGA) , 
several recent works highlight the importance of disorder, 
both positional and electronic. Stability of the coUinear 
ferromagnetic state has been investigated with randomly 
distributed Mn ions, and noncollinear ordering is sug- 
gested to be common to these semiconductor systems. 
[19] Competing (AF) interactions leading to frustration 
has already been evidenced by spin-glass behavior in II- 
VI DMS. [20] The presence of large compensation due 
to As antisite defects implies substantial electronic dis- 
order as well, and the sensitivity of Tc, magnetization 
M, transport and spin- wave spectrum to disorder has 
been investigated. [21-25] Monte Garlo simulations have 
also been used to study disorder effects on magnetic or- 
dering, [26-28] and dynamical and transport properties; 
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[29] the background fcrmions determine the spin interac- 
tions and hence the nature of the spin ordering, which 
in turn affects the fermionic states. Ab-initio methods 
[30-34] have also been recently employed. 

An alternative mechanism for the ferromagnetic cou- 
pling between impurity spins involves the hole-mediated 
RKKY interaction. [3] The RKKY theory has been ex- 
tended for various dimensionality structures, including 
effect of potential scattering through carrier mean-free 
path, indicating enhancement of ferromagnetic interac- 
tion by disorder in low dimensions. [7] Exchange and 
correlation has been shown to slightly enhance Tq within 
the RKKY theory. [9] Spin-wave dispersion in the RKKY 
picture has been compared with the result of spin-wave 
theory in which a uniform impurity-induced polarization 
has been assumed (VGA), resulting in a Zeeman splitting 
A in the carrier bands. [15] It was shown that the RKKY- 
level dispersion is incorrect except when A << Ep. 

The traditional RKKY approach is based on linear re- 
sponse in the weak-field limit (J << ep), which is not 
quite valid for the DMS. In this paper, we present a gen- 
eralized RKKY approach which takes into account the 
spatial variation of the impurity-induced carrier spin po- 
larization beyond linear response theory [section II]. In 
the generalized RKKY picture, the local magnetic field 
Bj = JSj of a magnetic impurity at site j polarizes the 
electrons locally, and the mobile band electrons spread 
this magnetic polarization in a characteristic manner: 
"^i = Xij{B)Bj, where Xij{B) represents the general- 
ized magnetic response. The spin Si of another mag- 
netic impurity placed at site i couples to this local elec- 
tronic magnetization, resulting in an effective generalized 
RKKY spin coupling J^Xij{J)Si-Sj. 

We find several interesting competing processes which 
limit the growth of spin couplings. As the RKKY re- 
sponse involves a particle-hole process, it vanishes for 
a filled (valence) band and grows with increasing hole 
concentration p. While the spin coupling is therefore ex- 
pected to strengthen with p, a competing process involv- 
ing the length scale sets in, which limits the growth of the 
spin coupling and therefore of the ferromagnetic transi- 
tion temperature T^. The Fermi wavelength Ap = ^w/k-F, 
which sets the RKKY oscillation length scale, decreases 
with hole doping, and therefore the spin coupling be- 
tween two magnetic impurities at a fixed separation goes 
through a maximum as a function of hole concentration 
[section III]. 

We find a similar optimization in the spin coupling as 
a function of the impurity field strength B. By going 
beyond linear response theory, and examining the gener- 
alized RKKY response for a fixed hole concentration, we 
find that the RKKY oscillation becomes more rapid with 
increasing polarizing field. Therefore, for a fixed separa- 
tion between two impurity spins, the spin coupling ini- 
tially increases like as expected, but then crosses over 
and eventually changes sign, resulting in frustration [sec- 
tion III]. The non- linear magnetic response thus brings 



out another limitation in the ferromagnetic spin coupling. 

In order to determine the extent to which the mag- 
netization response of a single impurity determines the 
macroscopic magnetic properties of the DMS, we have 
also considered a finite concentration of magnetic impu- 
rities distributed on a finite-size lattice [section IV]. Us- 
ing a novel Hubbard-JJ representation for the magnetic 
impurities in a DMS, we have studied the collective mag- 
netic response in the ferromagnetic state within a mcan- 
field-plus-spin-fiuctuation (MF-I-SF) approach. Treating 
the disorder aspects of the Mn-impurity system exactly, 
and electron correlation effects within the random phase 
approximation (RPA), our numerical analysis yields the 
spin stiffness from the low-lying collective (spin- wave) ex- 
citations [section VI] , which have a fundamental bearing 
on the ferromagnetic transition temperature Tc- Our ap- 
proach also allows for a quantitative study of the stabil- 
ity/instability of the Hartrcc-Fock (mean-field) ferromag- 
netic state [section V] , highlighting the presence of com- 
peting interactions. The Anderson Hamiltonian, with a 
hybridization term Vp^ between band fermions and the 
magnetic impurity orbital, has also been recently studied 
to obtain the ferromagnetic coupling between two mag- 
netic impurities. [35] 

II. MAGNETIC IMPURITY IN A HOST 

We consider a single-band spin-fermion lattice model 
-ff = ^ ekdk^^ak.a - ^ JiSi.&i , (1) 

k,(7 i 

with a double-exchange interaction between the magnetic 
impurity spin Si and the electron spin cTj at the impurity 
site i. The host (valence band) dispersion ck is taken 
to be parabolic for small k (top of the band at fc = 0), 
the k^ coefficient determining the inverse carrier mass 
m*. As the added holes go in long- wavelength states, 
the small-A; particle-hole processes near the Fermi energy 
are dominant, and therefore other details of the energy 
band are expected to be relatively unimportant. 

A. Host Green's function 

We consider an isotropic energy-band dispersion 
W 

efe = — cos ka (2) 

in three dimensions, with the wavevector magnitude ex- 
tending upto Tf/a. This dispersion incorporates the de- 
sired features, and yields a finite bandwidth without in- 
troducing sharp cutoffs. We choose length and energy 
units such that the lattice spacing a = 1 and the band- 
width W = 1. The advanced Green's function for the 
host is obtained as 
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FIG. 1. Real and imaginary parts of the local host Green's 
function. 
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Here 0{k) is a A;-space density of states, and for simplicity 
we choose a symmetric form 



= ak^ - bk^ 

= a{k - irf - b{k - 



(0 < fc < 7r/2) 

(7r/2 < fc < tt) , (4) 



so that the usual three-dimensional fc^ form is recovered 
for states near both the lower and upper band edges at 
k = IT and k = 0, respectively. We choose b = 2a /n'^, so 
that (}{k) is smooth at fc = 7r/2 (the slope dj3/dk = 0), 
and an overal normalization a = 120/77r^ so that the sum 
over states in the band J^^ 0{k)dk = 1. 

The above choice yields a symmetric band with a 
nearly semi-elliptical density of states, as seen in Fig. 
1, showing the real and imaginary parts of the local host 
Green's function go{u>)- Near the band edges, the real- 
part magnitude has a finite maximum and the imaginary- 
part has a square-root behaviour, as expected for the 
three-dimensional system. The band filling is shown in 
Fig. 2 as a function of the Fermi energy. 
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FIG. 2. Fermi-energy dependence of the band filling. 



l + (7Bjgo{uj) 



9ji{^) , (5) 



where go = gjj is the local host Green's function. The 
resulting local magnetization rrii at site i is then obtained 
as 



^Im[Gl(a;)-4(a;)], 

-no TT 



where 

Glico) - gU^) = gijiu) ATj gji{u;) , 
in terms of the T-matrix difference 

-2B, 



ATj = Tj - 



(6) 



(7) 



(8) 



Defining a field-dependent generalized magnetic re- 
sponse function Xij (B) through the relation 



m = Xij{B)Bj , 
Eqs. (6), (7), (8) yield 
du). 



(9) 



(10) 



B. Magnetic response 

We consider the impurity spin in the classical limit 
{JiSi Ji{Si) = Biz), and examine the magnetic re- 
sponse of electrons in a nearly filled band due to the 
magnetic coupling — ^jCTj.Sj, for an arbitrary strength 

of the impurity- induced local magnetic field Sj. 

For a single magnetic impurity at site j, the electronic 
Green's function G is exactly obtained in terms of the 
host Green's function g as 



III. EFFECTIVE SPIN COUPLINGS 

Another impurity spin Si placed at site i will couple 
with the local magnetization rrii produced by the local 
field of the spin Sj at site j, resulting in an effective 
interaction between the two spins given by 

ifspin( J) = -JiJjXij {B = JS)Si.Sj . (11) 
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FIG. 3. Magnetic susceptibility x(f ) for different hole dop- 
ing concentrations. 
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FIG. 4. Hole doping dependence of magnetic susceptibility 
for different r (corresponding to average Mn-Mn sepa- 
rations in a cubic host lattice with 5%, 8%, and 12.5% Mn 
impurity concentration) . 



A. Weak-coupling limit: RKKY interaction 

When the term in Eq. (10) can be neglected (vahd 
for B << W), one obtains a linear response 



rrii = XtjBj , 
where the magnetic susceptibility Xij 
duo 



(12) 



Xij 



\m.[gij{u;)gji{w)] 



= 4E E 



Em - 



(13) 



yields the standard oscillating RKKY interaction 

fl^RKKY = —JiJjXijSi-Sj . (14) 

The behaviour of Xij j sls a function of the separation r 
betwen the two sites i and j, shows that the oscillation 
becomes more rapid with doping (Fig. 3), as expected 
from the decreasing Fermi wavelength Ap = ^w/kp. 



Qualitatively similar results were obtained for a parabolic 
energy dispersion ek ~ fc^ with a finite bandwidth cutoff. 

For a fixed separation r/a = (l/x)^/"^, corresponding 
to the average Mn-Mn distance in a cubic lattice with 
Mn concentration x, the behaviour of xC?*) is shown in 
Fig. 4 as a function of band filling. The ferromagnetic 
coupling peaks at hole concentrations about 0.6 times the 
Mn impurity concentration. 



B. Generalized magnetic response 

It appears that the conventional RKKY picture based 
on the weak-coupling limit (B « W) cannot provide 
a good description of the interaction between Mn im- 
purities in Gai_2;Mn2:As. Core-level photocmission [36] 
yields J ^ 1 eV, which is comparable to the host band- 
width of « 2 eV for the heavy hole band. [37] It 
is therefore essential to go beyond the linear-response 
regime, and ior B ^ W we find that there are additional 
contributions in the generalized magnetic response func- 
tion Xij (B), which qualitatively modify the nature of the 
magnetic response and spin couplings. 



1. Impurity-state contribution 

For ui outside the band (|a;| > W/2), the T-matrix 

difference in Eq. (10) has imaginary terms of the type 
6{aj — Lu*), arising from the two poles 



l±Bgo{(^*)=0 



(15) 



corresponding to a spin-f impurity state at w| (below 
the lower band edge), and a spin-| impurity state at lo* 
(above the upper band edge). In three dimensions, 5o(i^) 
has a finite maximum at the band edges, and therefore 
impurity states are formed only when B exceeds a thresh- 
old strength B* . 

By expanding go(uj) near u>^, and expressing tJ as a 
simple pole, the impurity-induced correction is given by 



^ii ~ 9ii — 9ijT]gji — 



LO — LO^ — IT) 

where the impurity-state wavefunction tpl is given by 



(16) 



■sJ-dgo/du!\a,=u>^ 



(17) 



For any finite doping, only the spin-f impurity state is 
occupied, and the impurity-state contribution to the local 
magnetization is therefore simply obtained as 



(18) 



With increasing B, the impurity-state wavefunction be- 
comes more localized, and (ft — > 6ij as B — > oo. 
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FIG. 5. The r-dcpcndcnce of the generahzed magnetic re- 
sponse x(^i-B), for different field strengths. 




FIG. 6. The field-strength dependence of the generalized 
magnetic response x(^, .S), for a fixed hole concentration and 
Mn-Mn distance. 



2. Band contribution 

The other contributions to the imaginary part in Eq. 
(10) are from within the band (|a;| < W/2), and involve 
the real (imaginary) part of AT,- (u) and imaginary (real) 
part of gij{Lj)gji{uj). 

Including both the band and impurity contributions, 
the generalized magnetic response %('', B) evaluated from 
Eq. (10) is shown in Fig. 5 for different field strengths; 
the lowest-field case {B = 0.1) provides the RKKY re- 
sponse (B 0), for comparison. The length scale at 
which the first crossover from ferromagnetic to antifer- 
romagnetic coupling takes place is seen to decrease with 
increasing B. Fig. 6 shows the generalized magnetic re- 
sponse x{r^ B), for a fixed hole concentration and Mn-Mn 
distance. For small B = JS, the response is essentially 
constant (linear response), and in this regime the gen- 
eralized RKKY interaction energy i?^x(r, B) grows like 
J^, as in the mean- field and conventional RKKY pic- 
tures. However, the sharp suppression in the generalized 
magnetic response for B/W > 0.2 limits this growth and 
leads to a peak, which is seen to shift to higher B values 
with decreasing Mn-Mn separation (Fig. 7). This effect 
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FIG. 7. The generalized RKKY interaction energy 
B^xifj B) as a function of the local field strength B, for dif- 
ferent Mn-Mn distances. The host bandwidth is nominally 
taken as 10 eV. 



significantly increases the spin coupling in a higher Mn 
concentration system (such as GaMnN) beyond the fac- 
tor expected from the spin response (Fig. 4). 

As the effective carrier mass m* scales like the inverse 
bandwidth 1/W, the m* dependence of the generalized 
magnetic response x(r, B) can be directly deduced from 
Fig. 6, showing the B/W cx m* dependence for a fixed 
B — JS. The magnetic response in the (fixed) unit of 
1/B is obtained by multiplying xi^i B) in Fig. 6 (in unit 
of 1/M^) by B/W. This yields a linear m* dependence 
of x(r, B) (and hence the spin coupling energy and Tc) 
for low effective mass and then a sharp suppression with 
increasing m*. Sublinear dependence of at large m* 
has also been reported in Monte Carlo studies. [26] 

For a nominal host bandwidth of lOeV (with 1 eV 
PS 10^ K), the peak interaction energies are about 250 
K, 600 K and 1400 K for Mn concentrations of 5%, 8%, 
and 12.5%, and hole concentrations of 3%, 5%, and 8%, 
respectively [Fig. 7]. From this spin interaction en- 
ergy {JS)^Xi the ferromagnetic transition temperature 
Tc can be estimated within the spin-fluctuation theory. 
For a nearest-neighbour quantum Heisenberg model (in- 
teraction energy J) on a hypercubic lattice (coordina- 
tion number z), the transition temperature is given by 
Tc = T^^/fs!, somewhat lower than the mean-field value 
T^^F = js{S + l)z/3. [38] Here /,f = (l/iV)Ek(l - 
7k)~^ ^ 1 is a geometrical spin-fluctuation factor, where 
7k = (cos kx + cos ky + cos kz)/3 in three dimensions. 

As the effective RKKY interaction term between two 
spins is J'^Xij^i-Sj^ take J = J^x ^^id obtain Tc w 
2J'^xS{S + 1) for 2 = 6. Taking a realistic bandwidth 
oi W = 2eV for the heavy valence band, [37] the peak 
energy in Fig. 7 translates to a peak Tc of about 150 K 
and 850 K for 5% and 12.5% Mn concentrations, quite 
close to the observed highest Tc values for Gai_j;Mna;As 
and Gai_a;Mna;N. 
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IV. HUBBARD-(7 REPRESENTATION OF 
MAGNETIC IMPURITIES 



Wc now consider a (purely fermionic) Hubbard-model 
representation for the randomly distributed magnetic im- 
purities on a cubic lattice: 

H = t ^ (al^aja + h.c.) +t' ("L^jV + h.c. 



<lj>cr 



<IO>a- 



+ Cd ^ a\^aia + (n/T - n/) (n/^ - m) , (19) 

where I refers to the impurity sites, is the impurity on- 
site energy and nj = {ni-^ + nii)/2 is the spin-averaged 
impurity charge density. Higher spin magnetic impuri- 
ties, such as the S = 5/2 Mn impurities in Gai_xMnxAs, 
can be realistically represented within a generalized Hub- 
bard model representation involving multiple orbitals 
and different interaction processes (direct and exchange 
type, with respect to orbital indices). [39] For simplicity, 
we have taken the same; hopping (f/ = t = 1) bctwccm 
the host-host and host-impurity nearest-neighbour pairs 
of sites. The energy-scale origin is set so that the host on- 
site energy is zero, and wc take the impurity level to lie at 
the top of the host band (e^ = 6). The form of the Hub- 
bard interaction term is such that in the Hartree-Fock 
approximation it reduces to the double-exchange term. 



A. Hartree-Fock ferromagnetic state 

In the Hartree-Fock (mean-field) approximation, the 

interaction term reduces to a magnetic coupling of the 
electron to the local mean (magnetic) field A/: 
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FIG. 8. Instability of the undoped (HF) ferromagnetic 
state. The Goldstone mode ([/Ac = 1) corresponds to the 
minimum eigenvalue of the [x''(w = 0)] matrix, indicating 
maximal instability. 
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FIG. 9. Stabilization of the ferromagnetic state with hole 
doping; the Goldstone mode now corresponds to the maxi- 
mum eigenvalue. 

the MF Hamiltonian, and this procedure is iterated till 
self-consistency is achieved. 



where the electronic spin operator <?/ = 'I'j[(t]5'/ in terms 



of the spinor i 



ail 



, and the mean field A/ is self- 



consistently determined from the ground-state expecta- 
tion value: 



2A/ = U{ai) 



(21) 



Thus, in the classical (Hartree-Fock) limit, the inter- 
action term reduces to the corresponding form of the 
double-exchange term, with the mean field A/ represent- 
ing the impurity-induced local magnetic field Bj. 

Starting with an initial uniform mean field A/ = zA, 
the mean-field (MF) Hamiltonian is numcric;ally diago- 
nalized for a finite lattice to obtain the fermion eigen- 
functions (pia- and eigenvalues Ei^- The spin densities 
n/o- = X!^, <E-p^'i'iaf' yidd the new local mean fields 
Aj = C/(n/^ — nii)/2, which are then used to update 



B. Stability of the HF state 

The self-consistent, HF ferromagnetic state, with all 
local moments aligned in the same symmetry-breaking 
direction, does not necessarily represent a stable (lowest- 
energy) state. This is because the HF state really rep- 
resents an energy extremum, which may correspond to 
a saddle point having local energy minimum and max- 
imum along different directions in the order-parameter 
space. The stability of the HF state with respect to trans- 
verse perturbations in the order parameter is indicated 
by the maximum eigenvalue Amax of the [x°(ix' = 0)] ma- 
trix (Eq. 24). [40] The HF state is stable if C/Amax = 1, 
correspondng to the massless Goldstone mode, represent- 
ing a rigid rotation of the ordering direction. Instability 
is indicated if J/Amax > 1, signalling a growth of trans- 
verse perturbations about the HF state, which can also 
be interpreted as negative-energy bosonic modes. 
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Figure 8 shows some of the cigcnvahics of the [x''(cj)] 
matrix (inchiding the minimum and maximum) for the 
undoped HF ferromagnetic state of an 8"^ system, with a 
semi-ordered arrangement of 32 magnetic impurities (see 
section V for details). For ui = 0, the Goldstone mode 
{UXg = 1) is seen to correspond to the lowest eigen- 
value, indicating maximal instability of the ferromagnetic 
state. The structure of the eigenvector corresponding to 
the maximum eigenvalue indicates a tendency towards 
AF ordering of the impurity spins. The AF coupling 
arises from the exchange interaction J' ~ t'^/U due to 
the effective hopping t' (associated with impurity-band 
formation) between impurity sites. In the absence of the 
hole-induced (RKKY) ferromagnetic coupling, this ex- 
change interaction dominates and favours AF ordering 
of impurity spins. With hole doping, the ferromagnetic 
state gets stabilized, and the Goldstone mode now corre- 
sponds to the maximum eigenvalue (Fig. 9). 



C. Transverse spin fluctuations 

Transverse; spin fluctuations arc gapless, low-energy 
excitations in the broken-symmetry state of magnetic 
systems possessing continuous spin-rotational symmetry. 
Therefore, at low temperatures they play an important 
role in diverse macroscopic properties such as existence 
of long-range order, magnitude and temperature depen- 
dence of the order parameter, magnetic transition tem- 
peratures, spin correlations etc. 

We study the time-ordered, spin-wave propagator in- 
volving the spin- lowering {S^) and spin-raising {Sj~) op- 
erators at sites i and j: 

X..,it-t')^^{^o\T[Sr{t)S+it'mo) ■ (22) 

At the RPA level, the spin- wave propagator in frequency 
space is given by 



[x»] 



l-[U][x^{u,)] ' 



(23) 



where the zeroth-order, antiparallel-spin particle-hole 
propagator 

[x°H]iJ = i j ^Gl.{J)G\,{J - u:) (24) 



Bm>EF 



= E 



+ 



Ei<Ef 



Emi — Ell + ^ — Ell ~ ^ 



is evaluated using the eigenvalues Ei^^ and eigenvectors 
(/ijo- in the self-consistent, broken-symmetry state. In Eq. 
(23), the diagonal interaction matrix [U]ii = Udn has 
non-vanishing elements only at the magnetic impurity 
sites. For site-dependent interactions, it is convenient to 
recast Eq. (23) using simple matrix manipulations: 



where [A(a))] = [U] — [U][x^ {oj)][U] is a symmetric ma- 
trix, having non-vanishing matrix elements only in the 
reduced impurity basis: 



[Aico)]jj = U{l-U[x°i^)]u). 



(26) 



Spin-wave modes, represented by the poles in the prop- 
agator [x '"(w)], are hence given by the poles of the ma- 
trix [j4(a;)]/,7, as [U] is non-singular. In terms of the 
eigenvalues X„ and eigenvectors ipn of the [x°('-^)]/J 
trix, the spin-wave energies a;„ are therefore given by 

1 - UXniiOn) = . (27) 



V. SPIN- WAVE ENERGY 

The spin couplings and stiff'ness in the ferromagnetic 
state can be determined from the spin-wave energies. To 
see how the magnitude and sign of the spin couplings de- 
pend on impurity separation, we have considered several 
impurity arrangements with different numbers (iVimp) of 
magnetic impurities in a cubic host lattice with N = 8^ 
sites. These arrangements include: (i) an ordered im- 
purity arrangement of 64 impurities (x = 1/8) on a cu- 
bic superlattice with impurity separation 2a, (ii) a semi- 
ordered arrangement of 32 impurities (x « 6%), with 
same NN impurity separation (2a) in the z direction but 
greater in-plane separation (v^a), and (iii) a disordered 
arrangement of 30 impurities (x « 6%) with NN separa- 
tions ranging between 2a and 3a. 

For the undoped (insulating) state, we take = N 
and Ni = N—Ni^p-, all spin-| impurity states (pushed up 
by the local mean field) are then unoccupied, resulting 
in local-moment formation. Hole doping is introduced 
by reducing , and band fillings are so chosen that the 
Fermi energy lies in gaps between (nearly) degenerate 
groups of eigenvalues. 

The undoped self-consistent ferromagnetic state is 
found to be maximally unstable, as discussed earlier. In- 
deed, the self-consistent antifcrromagnetic state is actu- 
ally found to be stable, confirming the dominance of the 
AF spin couplings ,/' ~ t''^ /U. With hole doping, the fer- 
romagnetic state is stabilized, and the spin-wave energies 
LOn are extracted from the pole condition UXn{i^n) = 1- 
Hole doping dependence of the lowest spin-wave energy 
wiow is shown in Fig. 10 for the ordered (i) and dis- 
ordered (ii) impurity arrangements, with U = 4 and 5, 
respectively. Figure 11 shows the ?7-depcndence of the 
lowest spin-wave energy for the ordered (i) and disor- 
dered (iii) arrangements, with = 482 (p fa 6%) and 
N^i = 505 {p ~ 1.4%), respectively. The optimization of 
the spin coupling with respect to both hole doping and 
interaction energy U is qualitatively similar to that in 
the RKKY picture. 



[x-+(c^)] 



[A{u;)] [U] 



(25) 
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FIG. 10. The lowest spin-wave energies for the ordered (+) 
and disordered (x) arrangements peak near 60% fractional 
hole concentration, very similar to the RKKY spin coupling 
J^Xij (Fig- 4). Increasing impurity separation lowers the spin 
stiffness. 
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FIG. 11. Optimization of the lowest spin-wave energy with 
U , qualitatively similar to the generalized RKKY spin cou- 
pling (Fig. 7). 

While the lowest spin-wave energy is softened by dis- 
order, the highest spin-wave energy is, however, signifi- 
cantly enhanced, as shown in Figure 12. This enhance- 
ment of Whigh is associated with localization of spin- wave 
states over impurity clusters in which the relatively closer 
spins are more strongly coupled. [44] For the same min- 
imum impurity separation (2a) in arrangements (i) and 
(ii), disorder-induced localization leads to stronger bonds 
between the cluster spins. Also shown (for the ordered 
case) is the Stoner (single-particle excitation) gap, which 
is roughly proportional to the MF impurity magnetiza- 
tion. The spin-wave branch merges with Stoner excita- 
tions at about 8% hole concentration. 

For the ordered impurity arrangement, the spin- wave 
energy range allows the spin couplings to be extracted, as 
discussed below. Assuming nearest-neighbour exchange 
interaction J between the impurity spins on the super- 
lattice, the spin- wave energies are given by 



Wq = JSz{l - 7q) , 



(28) 



FIG. 12. The highest spin-wave energy for the ordered (-I-) 
and disordered (x) impurity arrangements, showing disor- 
der-induced stiffening of the high-energy mode. 

where 7q = (cosq^ -|- cos qy + cos qz)/'S. The spin- wave 
modes on the impurity super lattice are plane waves, with 
wave- vector components given by = n^27r/i, where 
are integers and L = 4 for the 64-impurity super- 
lattice. The wave vectors q ~ (l,0,0)27r/L etc. and 
q = (2, 2, 2)2n/L correspond to the lowest- and highest- 
energy modes, respectively. The corresponding ener- 
gies wiow = JSz/3 and t^high = 2JSz yield a ratio 
Whigh/t^iow = 6. We not only find the actual ratio to 
be quite close (about 7 for most doping cases), but the 
degeneracies in the spin-wave spectrum are also in close 
agreement, indicating that nearest-neighbour spin cou- 
pling is dominant. 



VI. CONCLUSIONS 

A comparitive study of a generalized RKKY approach 
and a MF-I-SF approach offers new and useful insight 
into the mechanism of carrier-mediated ferromagnetic or- 
dering in a dilute magnetic semiconductor. While the 
MF+SF approach provides quantitative understanding 
of the spin couplings, competing interactions, spin- wave 
excitations, low-temperature spin dynamics, and the crit- 
ical temperature, the generalized RKKY approach pro- 
vides a qualitative understanding in terms of a sim- 
ple physical picture involving the impurity- induced oscil- 
lating carrier-spin polarization, which complements the 
MF-I-SF approach. Our key finding is an optimization of 
the spin coupling (spin- wave energy) with respect to both 
hole doping and the impurity polarizing field strength (J 
or U), which is in agreement with recent Monte Carlo 
study, [28] and can be physically understood in terms of 
a competition between the increasing magnitude of the 
carrier-spin polarization and increasing rapidity of its os- 
cillation. We find that the optimum (fractional) hole con- 
centration for the spin coupling occurs at p/x « 0.6, and 
both the spin coupling energy J'^x{r, J) and the spin- 
wave energy scale with the carrier bandwidth W, for 
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fixed J/W or U/W. The oscillating spin polarization 
also highlights the role of competing interactions in the 
instability of the collinear ferromagnetic state. 

In this paper, we have presented the first study of spin- 
wave excitations in the ferromagnetic state of the DMS 
within a microscopic correlated lattice fermion model 
which treats finite impurity concentration, impurity dis- 
order, and electron correlation on an equal footing. With 
regard to electron correlation, the MF+SF approach has 
been shown to be applicable in the full range from weak 
to strong coupling, [41] and extensively used in the con- 
text of strongly correlated layered cuprate antifcrromag- 
nets which exhibit pronounced spin fluctuations. When 
the spin- wave energy is much smaller than the mean-field 
strength A, spin dynamics is dominant at low tempera- 
tures and charge fluctuations can be ignored for T << A. 
Incorporating the low-energy spin fluctuations about the 
MF state yields quantitatively correct temperature de- 
pendence of (sublattice) magnetization and reliable Tq 
within the rcnormalizcd spin- wave-theory (SWT). [41] 
Whereas Tq pertains to global ordering, with the spin 
coupling energy providing the relevant energy scale for 
spin fluctuations, the MFT deals with local ordering, 
and greatly over-estimates the transition temperature 
(Tc ^ A), which really represents the moment-melting 
temperature. 

Specifically with regard to the DMS, there is a subtle 
issue concerning the energy scale relevant for global or- 
dering. Whereas for a generic ferromagnet, energy scales 
corresponding to the local mean field and spin coupling 
arc identical, for the DMS, three distinct energies can be 
identified — the two local mean fields seen by the car- 
rier spin (~ J) and impurity spin (~ J^Xii)) and the 
coupling between impurity spins J'^Xij)- In the weak 
doping limit {p/x 0), the magnetic response function 
X{r) decays slowly on the impurity-separation scale, so 
that Xij ~ Xii) and the distinction between the two latter 
energy scales gets blurred. However, for a realistic frac- 
tional doping of p/x ~ 15%, the impurity spin couphng 
J^Xij is by far the lowest energy scale, and should there- 
fore control the low-temperature behaviour of the mag- 
netization M{T). The ferromagnetic transition temper- 
ature Tc, determined within the spin- wave- theory from 
the spin coupling energy for a realistic (heavy) hole band- 
width, corresponds closely with the observed Tc values in 
Gai-xMuxAs and Gai_xMnxN. With appropriate hole 
doping, transition temperature much above room tem- 
perature appears possible for x = 1/8, which is within 
experimental limit. [42] 

The MF-I-SF approach also highlights the role of impu- 
rity disorder. While the low-energy spin-wave modes are 
significantly softened as compared to the ordered case, 
the high-energy spin-wave modes arc clearly stiffened, 
indicating that a single spin-wave energy scale is not suf- 
ficient to describe the low-temperature spin dynamics in 
the DMS. In fact, a distribution in spin couplings, with 
weak and strong bonds, has been suggested to be re- 
sponsible for anomalous temperature dependence of mag- 



netization, susceptibility, specific heat etc. [43] Using a 
simple model involving two spin-excitation energy scales 
corresponding to weakly and strongly coupled spins, the 
temperature dependence of magnetization is found to be 
in good agreement with the SQUID magnetization data 
for Gai_a;Mnj,As. [44] 

ACKNOWLEDGEMENT 

Helpful discussions with T. Pareek and R. C. Budhani 
are gratefully acknowledged. 



^ H. Ohno, H. Munekata, T. Penny, S. von Molnar, and L. 

L. Chang, Phys. Rev. Lett. 68, 2664 (1992). 
^ H. Ohno, A. Shen, F. Matsukara, A. Oiwa, A. Endo, S. 

Katsumoto, and Y. lye, Appl. Phys. Lett. 69, 363 (1996). 
^ F. Matsukara, H. Ohno, A. Shen, and Y. Sugawara, Phys. 

Rev. B 57, R2037 (1998). 
" N. Theodoropolpu et al., AppL Phys. Lett. 78, 3475 

(2001) . 

5 M. L. Reed et al, AppL Phys. Lett. 79, 3473 (2001). 

S. Sonoda et al, J. Cryst. Growth 237-239, 1358 (2002). 

T. Dietl, A. Haury, and Y. M. d'Aubigne, Phys. Rev. B 

55, R3347 (1997). 
* M. Takahashi, Phys. Rev. B 56, 7389 (1997). 
^ T. Jungworth, W. A. Atkinson, B. H. Lee, and A. H. 

MacDonald, Phys. Rev. B 59, 9818 (1999); B. H. Lee, T. 

Jungworth, and A. H. MacDonald, ibid 61, 15606 (2000). 
^° T. Dietl, H. Ohno, F. Matsukara, J. Cibert, and D. Fer- 

rand. Science, 287, 1019 (2000). 
" T. Dietl, H. Ohno, and F. Matsukara, cond-mat/007190 

(2000) . 

T. Dietl, F. Matsukara, and H. Ohno, Phys. Rev. B 66, 
033203 (2002). 

T. Jungworth, J. Konig, J. Sinova, K. Kucera, and A. H. 

MacDonald, Phys. Rev. B 66, 012402 (2002). 
" H. Ohno, N. Akiba, F. Matsukara, A. Shen, K. Ohtani, 

and Y. Ohno, Appl. Phys. Lett. 73, 363 (1998). 

J. Konig, H. H. Lin, and A. H. MacDonald, Phys. Rev. 

Lett. 84, 5628 (2000); J. Konig, T. Jungworth, and A. 

H. MacDonald, Phys. Rev. B 64, 184423 (2001). 
^® A. Chattopadhyay, S. Das Sarma, and A. J. Millis, Phys. 

Rev. Lett. 87, 227202 (2001). 
^'^ S. Das Sarma, E. H. Hwang, and A. Kaminski, cond- 

mat/02 11496 (2002). 

G. Bouzerar and T. P. Pareek, Phys. Rev. B 65, 153203 

(2002) . 

^® J. Schliemann and A. H. MacDonald, Phys. Rev. Lett. 
88, 137201 (2002). 

^° J. K. Furdyna and J. Kossut, Semiconductors and 
Semimetals, Vol. 25 (Academic Press, New York, 1988). 
21 M. Borciu and R. N. Bhatt, Phys. Rev. Lett. 87, 107203 

(2001) . 



9 



^■^ A. L. Chudnovskiy and D. Pfannkuche, Phys. Rev. B 65, 
165216 (2002). 

23 M. Berciu and R. N. Bhatt, Phys. Rev. B 66, 085207 
(2002). 

2* G. Bouzerar, J. Kudrnovsky, and P. Bruno, cond- 

mat/0208596 (2002). 
2^ C. Timm, F. Schafer, and F. von Oppen, Phys. Rev. 

Lett. 89, 137201 (2002); C. Timm and F. von Oppen, J. 

Supercond. 16 23 (2003); cond-mat/0209055 (2002). 
^® J. Schliemann, J. Konig, and A. H. MacDonald, Phys. 

Rev. B 64, 165201 (2001). 
2'' M. P. Kennett, M. Berciu and R. N. Bhatt, Phys. Rev. 

B 66, 045207 (2002). 
2* G. Alvarez, M. Mayr, and E. Dagotto, Phys. Rev. Lett. 

89, 277202 (2002). 
2" G. Alvarez and E. Dagotto, cond-mat/0303350 (2003). 

H. Akai, Phys. Rev. Lett. 81, 3002 (1998). 
^1 M. Shirai et al, J. Magn. Magn. Mater. 177-181, 1383 

(1998). 

32 S. Sanvito, P. Ordejon, and N. A. Hill, Phys. Rev. B 63, 
165206 (2001). 

33 L. M. Sandratskii and P. Bruno, Phys. Rev. B 66, 134435 
(2002). 

3* L. Craco, M. S. Laad, and E. Miiller-Hartmann, cond- 

mat/030342 (2003). 
3^ V. A. Ivanov, P. M. Krstajic, F. M. Pootors, V. Fleurov, 

and K. Kikoin, Euro. Phys. Lett. 61, 235 (2003). 
3® J. Okabayashi et al, Phys. Rev. B 58, R4211 (1998). 
3^^ J. R. Chclikowski and M. L. Cohen, Phys. Rev. B 14, 

556 (1976); W. A. Harrison, Electronic Structure and the 

Properties of Solids (Freeman, San Francisco, 1980). 
3« M. E. Lines, Phys. Rev. 135, A1336 (1964); S. Doniach 

and E. H. Sondheimer, Green's Functions for Solid State 

Physicists (Benjamin, New York, 1974). 
3» A. Singh and P. Sen, Phys. Rev. B 57, 10598 (1998). 

For a stability analysis of the hole and electron doped 

antiferromagnetic state, see A. Singh and Z. Tesamovic, 

Phys. Rev. B 41, 614 (1990); A. Singh and H. Ghosh, 

Phys. Rev. B 65, 134414 (2002). 
*^ A. Singh, cond-mat/9802047 (1998); Euro. Phys. Jour. B 

11, 5 (1999); and references therein. 
*2 A. ,L Blattner and B. W. Wessels, cond-mat/0205602 

(2002). 

■'3 M. P. Kennett, M. Berciu and R. N. Bhatt, Phys. Rev. 
B 65, 115308 (2002). 
A. Singh, cond-mat/0307009 (2003). 



10 



